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We study holographic thermalization in spacetimes with a chemical potential and a non¬ 
trivial dilaton field. Three non-local observables are used to probe the whole process and 
investigate the effect of the ratio of the chemical potential over temperature x and the 
dilaton-Maxwell coupling constant a. It is found that the saturation time is not always 
a monotonically increasing function of y, the situation depends on a. When 0 < a < 1, 
larger % yields longer saturation time, while for a > 1, the situation becomes more complex. 
More interesting, we found that although a indeed has influence on the whole thermalization 
process, it nearly does not affect the saturation time, which indicates the universality of the 
saturation time for the dual one-parameter field theories. 
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I. INTRODUCTION 


In the past two decades, one of the most important discoveries in modern theoretical physics 
is the AdS/CFT correspondence relating strongly coupled field theories to weakly coupled 

gravity. It provides a powerful tool to study strongly coupled field systems where traditional 
methods encounter great challenges or even break down. Fruitful results have been obtained bv 
applying this duality to various areas of modern physics, such as condensed matter physics jd-fj 
and QCD 13H2|. 

One noticeable application of this duality is to study non-equilibrium physics, in particular 


in the so-called ’’holographic thermalization” topic 


10-[45|] which is an aternpt to describe the 


thermalization process of the quark-gluon plasma (QGP) produced in heavy ion collisions at the 
Relativistic Heavy Ion Collider (RHIC) {§, 46-48]. Although the near-equilibrium physics of the 


QGP can be understood using a hydrodynamic description 


49-1521]. the whole process to reach 


thermal equilibrium is yet not known. Observed data from RHIC show that the time scale needed 
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for QGP to reach thermal equilibrium is considerably shorter than expected based on perturbative 


methods 


3 , 54]. 


This indicates that the thermalization process of QGP is strongly coupled, and 


thus motivates us to use AdS/CFT relation. 

According to AdS/CFT, it is proposed that this thermalization process is dual to the process 
of a black hole formation in the bulk via gravitational collapse of a scalar field in AdS space [10- 


[jjj], or t° a collapse of a thin-shell matter described by an AdS-Vaidya metric [I5l4l8l]. To probe 
the detailed process of the thermalization, local-observables in the boundary field theory, such as 


energy-momentum tensor and its derivatives, are not sufficient. In refs. 


17 


181 ]. in the model 


of the collapsing of a thin-shell, three non-local observables have been proposed to probe the 
process, the two-point function, the Wilson loop and the entanglement entropy. According to the 
AdS/CFT relation, in the semiclassical limit these three non-local observables in the boundary 
field theory can be evaluated by calculating certain geometric quantities in the bulk geometry. 
By studying the time evolution of these three non-local observables, and under Einstein gravity 
frame, the thermalization process has thus been studied in detail and some interesting results are 
found. First, the saturation time which is needed to reach thermal equilibrium depends on the 
geometric size of the probe in the boundary field theory. Probes with smaller size thermalize 
faster indicating a ’’top-down” thermalization mechanism, which is very much in contrast with 
the standard ”bottom-up” paradigm [531 ] from the use of perturbative methods. Moreover, all 


probes used show a delay in the onset of the thermalization, an apparent non-analyticity in the 
end of thermalization, the transition to full thermal equilibrium is instantaneous and these features 
are independent of dimensionality. Finally, for homogeneous initial conditions the entanglement 
entropy thermalizes slowest, and thus sets a timescale for equilibration that saturates a causality 
bound. 

Later, this study is extended to include the effect of chemical potential in the boundary field 
theory 23], |24J (see also 33y), which is usually the case in real heavy ion collision processes. In 
the gravity side, this is dual to introducing a gauge field and considering charged background 
geometries. It is found that the time evolution of the three non-local observables exhibit similar 
behavior as in chargeless case. The arise of the chemical potential affects the saturation time: the 
larger the chemical potential, the longer the saturation time. Also the study has been extended to 
other situations by considering hig her curvature corrections ([27, 31, 32], non-linear electromagnetic 


effect 


441 ], angular momentum [30j], noncommutative 4l| and hyperscaling violating geometries 


ies Q, 


361 ]. Further, it is extended to study holographic thermalization of field theories living in curved 


spacetimes [34 


, 37, 3, 45]. 
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On the other hand, it is believed that Einstein gravity should be considered as an effective 
description of the underlying quantum gravity theory (string theory for example) at low energy. 
Thus, it is interesting to study the effect of stringy corrections on the thermalization process. 
As we have mentioned above, some work has been done in this directions by considering higher 
curvature corrections 271.131. l32l. 1451]. It is found that the existence of higher curvature corrections 
considerably affects the thermalization i process. For example, by including the effect of Gauss- 


Bonnet terms, the authors in refs. [271 . |31 


451 ] found that a larger Gauss-Bonnet coupling 


constant will make the saturation shorter, thus leaving the imprint of this stringy effect. In this 
paper, we would like to consider another kind of stringy effect, the emergence of a dilaton field. The 
action we consider is the well-known Einstein-Maxwell-dilaton action which is a low-energy effective 
action of string theory [53], where the dilaton couples with the Maxwell field in a non-trivial way. 


To use the AdS/CFT relation and study the effect of the dilaton field, we should first find asymp¬ 
totically AdS solutions with non-trivial dilaton configuration under this frame. It is a non-trivial 
task. Amounts of work have been done to explore various solutions of this gravity theory. Non¬ 


trivial coup 
spacetimes 


ing between the dilaton and the gauge field yields interesting effects on the available 


56-1661]. However, all of these black hole solutions with non-trivial scalar configuration 


are not aysmptotically AdS. In [67 69(], the authors even argued that, with only one Liouville-type 


lO 

e so 

70 

m 


dilaton potential besides the cosmological constant, there are no dilaton-dS or AdS black 
lutions. After lots of efforts, an aysmptotically AdS black hole solution was found in ref. 

By choosing an appropriate combination of three Liouville-type dilaton potentials, the authors 
successfully obtained the static dilaton black hole solutions which are asymptotically (A)dS in 


72h 75I] and has shown 


four and higher dimensions. This solution has some interesting properties 

potential applications in cosmology 7 ( 3 ]. We will work with this charged dilaton AdS black hole, 
and see if and how the presence of the non-trivial dilaton field affects the thermalization process. 

The paper is organized as follows. In the next section, we briefly introduce the charged dilaton 
AdS black hole and then generalize it to obtain its Vaidya-like version. In Sec. Ill, by applying the 
three non-local observables, we study the detailed process of thermalization to see the influence of 
the chemical potential and the dilaton field. The last section is devoted to summary and discussions. 

II. BULK SPACETIMES 


In this section, we briefly introduce the charged dilation black hole in AdS space, and obtain 
its Vaidya-like version by introducing external matter source. 




























4 


A. Charged Dilaton Black Hole in Anti-de Sitter Space 

We consider the action of (n + l)-dimensional Einstein-Maxwell-dilaton gravity 0 , fl, Ed, E 3 

5 = “ lL J r+lx ^ ( R ~ , ( 1 ) 

where a is the coupling constant between the dilaton and the Maxwell held. When a = 0 this 
action reduces to the usual Einstein-Maxwell-scalar theory, while for a = 1 the last term in the 
bracket represents the well-known dilaton-Maxwell coupling that appears in the low energy string 

59]. The dilaton potential is chosen to take the form 


action in Einstein’s frame 
2A 




, . ( - a 2 [(n + l) 2 - (n + l)a 2 - 6(n + 1) + a 2 + 9] e -4(n-2)®/[(n-i)a] 

+(n-2) 2 (n-a 2 )e 4Q<3i /( ri - 1 ) +4a 2 (n- l)(n-2)e- 2 ( n - 2 - Q2 ^/[( n - 1 ) Q ]|, (2) 

where A = — - 7 ^ 2 —- is the cosmological constant and l the AdS radius. This type of dilaton potential 

can appear in the compactihcation of a higher-dimensional theory, including various supergravity 

models in four dimensions [t 7 , 3- From the above action, we can see that reversing the sign of a 

is equivalent to reverse the sign of the dilaton. Thus, without loss of generality, we can restrict a 

to be non-negative, a > 0 . 

The corresponding equations of motion are 

1 4 

R ' w ~ 2 9 ^ R ~ 


- l -9„v (V<f>) 2 - 

- 2e -4«*/(n-i) (f, ip FJ‘ - jg^FpaFlA = 0, 
V 2 ^ n-\ d\ a 4 a ^,/( n _i) pa _ 

v 8 a$ + 2 pa ’ 

V M ^ e — 4 «*/(»—1 )f^ = 0. 


( 3 ) 


They admit a charged dilaton black brane solution 

dp 2 


3, Ei] 


ds 2 = -N(p)f(p)dt 2 + jj-j+p 2 g(p)dx 2 n _ 1 , 


N ( P ) = T-^- 3 \ f( P ) = ^[n-2)^ __ __T 1-7 , 


g{p) = t 7 


P 

T = 1 - 


71 — 2 


( 4 ) 


where the constant 7 = ( n „ 2 )(?T-?+Q a ) • The dilaton and the electromagnetic fields take the form 


n—1 , _ _ e 4a$/(n-l) 

$ ( r ) = 4 VtI 2 + 27 - "7) ln T, 


( 5 ) 
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The remaining components of vanish. The parameters m and q are related to the physical 
mass M and charge Q of the black brane, 


M = —t~~ 1 (n — l)m, 


Q = — ;— Q, 


( 6 ) 


167T V ' ' ^ 4 TT 

where Q n -i is the volume of dx 2 _ l . Hereinafter, without risk of confusion, we simply call the 

parameters m and q the mass and charge parameter respectively. The constant b is related to the 

mass and the charge parameters through the relation 

2 (n — l)(n — 2) 2 _ 2 

q 2(n — 2 + a 2 ) m ‘ 

We should emphasize here that the solution contains three free parameters: ( a,m,q ), or equiva¬ 
lently ( a,m,b ). 

When a = 0, the solution reduces to the well-known Reissner-Nordstrom-AdS black brane. The 
holographic thermalization in this case has been thoroughly discussed in refs. J 23 , 


241 ]. Thus, in 


the following we will mainly focus on the case with a > 0. In this case, we note that the solution 
is not real for 0 < p < b. We should exclude this region from spacetime, what can be achieved by 
defining a new radial coordinate r as [tI 74 ] 


The metric becomes 

„2 


ds z = —N(r)f(r)dt z + 


2 2 ,2 

r = p - b . 


-r*2 

2 | 1 hi 1 / 2 , l2\ / \ j-»2 


( 8 ) 


(r 2 + b 2 )f(r) 


+ {r +b )g(r)dx n _ 1 , 


N(r) = r~( n ~ 3 ^, f( r ) = b2 r( n ~ 2 ^ - 


m 


g(r) = r\ 


V 2 

r = 1 - 


( r 2 + fr2)(n—2)/2 
n —2 


pl-7 


(9) 


. y/r 2 + b 2 

with the coordinate r now valued in the range 0 < r < 00 . The dilaton and the electromagnetic 
field become 

y/]V(r)re4“ $ /( n i) 


n — 1 


^( r ) = —:— V 7( 2 + 27 - nj) InT, 


Ffr — — - 


( 10 ) 


4 v ' v 1 1 L ' g(r)( n— -^)/ 2 (r2 -)- f>2) n /2 

It can be shown that r = 0 is a curvature singularity. The event horizon is the largest root of 
equation /(r^) = 0. Depending on the free parameters (a,m,b) (or equivalently (a,m,q)), there 
can be two horizons (inner and outer horizons), one degenerate horizon (extreme case) and naked 
singularity. The Hawking temperature on the event horizon is 

(JV/)' 


Th = x [F+ d 


4i t^/Nt 


( 11 ) 


r=r h 
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It is not possible to give an analytical expression of r^. However, we can express the mass parameter 
m in terms of 77 by solving f{rh) = 0 , which is 

( r h + b2 T /2 r (n—1)7—i 


m = 


l 2 


( 12 ) 


with Fh = r(r = rh). Substituting m into Eq. (jTTJ) , one obtains the temperature as 


T h = 


bf 


l 7 (n-l)/ 2 -l 

h 


-|(n- 2)[(n - 1)7-1] - (n- 1) [(n - 2)7 - 2]r h |. 


(13) 


4vr/ 2 (l - 

From it, we can see that when 1 > 7 > |(ct > 1 ), the temperature is always positive Tjj > 0 
and no extremal limit exists. While 0 < 7 < |(a < 1), like a Reissner-Nordstrom black hole, the 
temperature can be positive (non-extremal case) or zero (extremal case). In the latter case, if we 
think of all the parameters but b fixed, then the maximum value of b corresponding to the extremal 
case can be obtained by solving T#(&, rh,a ) = 0, which is 


t - =r 4( (n ~ 1)p r (n ~ 2)1 ) 



(14) 


For more discussions on the thermodynamical properties of this black brane, see refs. 


mm>\. 


According to AdS/CFT, the above black brane solution is dual to a CFT on the boundary with 
temperature given by Eq. ED- Moreover, the chemical potential /i of the boundary CFT is related 
to the asymptotical value of the temporal part of A^, i.e., n ~ lim^oo A t . By integrating Eq. ED, 
we can obtain the expression of At, 

^/{n — 1 )m&( n-2 )/ 2 


At = — 


+ Tr 


y/2 (n — 2 + a 2 )(r 2 + 6 2 )( n 2 )/ 2 
Where To is a constant corresponding to the electrostatic potential at r 
such that the gauge field vanishes at the horizon, i.e., 

•^/(n — 1 )mM n_2 )/ 2 

To = 


(15) 

00 , which is defined 

(16) 


yj2(n — 2 + a 2 )(r% + b 2 )( n 2 )/ 2 
The precise relation between the chemical potential and the gauge field is subtle because there is 
a dimension mismatch between the chemical potential in the dual field theory and the gauge field 
A^ defined in the action Eq. ©• Thus, one has to rescale the gauge field as A^ = where £ is a 
scale with length unit that depends on the particular compactification. At this time, the chemical 
potential and the gauge field has the same unit and the duality relation is 

To 

H = lim A t = —. 

r—>• 00 c 


( 17 ) 
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Since the boundary field theory we considered is conformal, the only relevant parameter we can 
vary is a dimensionless ratio x constructed from the chemical potential and the temperature, i.e., 
H 


t h ’ 


X ; 




d'n 


(18) 


T H T h £, 

Now we can see that if a and m are kept fixed, then by varying 6 which is equivalent to varying 
the charge parameter q, we can explore various values of the ratio x hi the dual field theory. For 
0 < a < 1, when b varies from 6 = 0 (vanishing To) to b ex t (vanishing Th), X spans the whole 
range of values, i.e., from y = 0 to oo. However, for a > 1, the situation becomes complicated. 
As we will see later, in this case, as 6 runs over the whole range of values, x is bounded in a finite 
range of values. 


B. Vaidya-like Solution 


Our goal is to investigate the thermalization process of the boundary field system under a 
quantum quench. According to AdS/CFT, this process can be holographically modeled by the 
process of a black hole formation in the bulk which can be described by a Vaidya-like metric. 
Thus, in this section, we would like to get the Vaiya-like metric by generalizing the black brane 
solution to be time-dependent. 

It is more convenient to work with a new radial coordinate £ = -y such that the AdS boundary 
lies at z = 0. Thus, the metric becomes 


ds 2 

N(z) 

9(z) 


7= 

.2 L 


-N{z)f{z)dt 2 + 
r _ (n _ 3 ) 7) /( 2 ) = 


dz 2 


z 4 (Z 4 + b 2 z 2 )f{z) 
Z 4 + b 2 z 2 


Z 4 + b 2 z 2 2 

H-o- g{z)dx n _ i, 




l 2 z 2 

tv r = i- 


r (™- 2 b _ m 
bz N n ~ 2 


\/Z 4 V b 2 z 2 


n —2 


r i~j 


Vl 4 + b 2 z 2 


(19) 


Introducing ingoing Eddington-Finkelstein coordinate 

Z 4 dz 


dv = dt — 


z 2 y/(l± + b 2 z 2 )N{z)f(z)' 


the above metric becomes 


( 20 ) 


ds 2 

f(z) 


—N(z)f(z)dv 2 - 


2 Z 4 


N(z) Z 4 + b 2 z 2 2 

dvdz H- k - g{z)aXn-i, 


Z 4 + b 2 z 2 
l 2 z 2 


p(n—2)7 _ 


m 


Z 4 + b 2 z 2 

z 

v/Z 4 + b 2 z 2 


n —2 


r 1-7 . 


( 21 ) 
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Usually, the Vaidya-like version can be obtained by using time dependent mass and charge param¬ 
eters (m —» m(v), q —>• q(v)), that is, 

m(v) = mT(v), q(v) = gT(u) 1 / 2 , (22) 


where m and q are the values of the corresponding parameters at late times. We can still introduce 
the parameter b, which satisfies the relation Eq. (ED, to make our expression explicitly. After taking 
such time-dependent parameters, we obtain the Vaiya-like metric which takes the same form as 
Eq. (1211) but with f(z) being now 


/4 i_ u2 / 

f(v,z) = —±__r(- 2 ^ - mT[v] 


VI 4 + b 2 z 2 


n —2 


r 1-7 . 


(23) 


Now, the equations of motion Eqs. ((3]) do not hold any more and external matter sources need to 
be introduced, 


R„u - \g^R - ^-j- (d^dv* - (V4>) 2 - ^-Vu(</>)) 

- 2e -4«*/(n-i) (f iip F u p - -^ hlu F fta F pa \ 


Svr T$*\ 
SnJexf 


(24) 


By checking the equation of motion of the dilaton field, we find that it always holds and no more 
additional external terms need to be introduced. 

In order to keep simplicity, in this paper we only consider the case with n = 4, that is, the 
dual field theory we considered is a four-dimensional conformal field theory. The above Vaidya-like 
metric is a solution of the equations of motion Eqs. 0) provided the external matter source satisfies 


8vr T (ext) = 

flU 


^ n J(ext) 


l 4 


3z 3 [(2 + a 2 )l 4 + a 2 b 2 z 2 
2(2 + a 2 )(l 4 + b 2 z 2 ) 5 / 2 \l 4 + b 2 z 2 

V 4 


3a z 

2(2 FF) 


mt{v) V 8 V 8\ 




V3bz 5 


\^V 2 ( 2 + « 2 ) V^ 4 + b 2 z 


2+ol z 


m- 



y/T(v) 

where dot means d v . We can see that the in-falling matter is charged dust. As in refs, 
consider T(v) to take the form 


(25) 


3,Q, 


we 


T(v) = - (1 + tanh— ], (26) 

2 V v oJ 

where vq is a small parameter, with typical order of 10“ 3 . That is, the ingoing falling matter is a 
thin shell. Thus, in the bulk gravity side, we consider the collapsing of an ingoing falling thin-shell 
to form a black hole. And this picture is dual to a sudden energy injection into the field system 
evolving towards thermal equilibrium. 
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III. NON-LOCAL OBSERVABLES 


Now, we will use the Vaidya-like metric constructed in last section to discuss the thermalization 
process of the dual field system. To investigate the details of the thermalization, one can use three 
non-local observables, the two-point coorelation function, the Wilson loop and the entanglement 
entropy, to probe the process. According to AdS/CFT, in the saddle approximation these three 
non-local observable can be evaluated by calculating the corresponding geometric quantities in the 
bulk geometry. 

As we stated above, we mainly focus on the case with n = 4. In doing numerical calculations 
in the following, we would like to set l = m = 1. That is, we fix the mass of the black brane 
formed finally, which is equivalent to fix the amount of energy injected into the dual boundary held 
system. And then we would like to see the effect of the chemical potential (more precisely y) and 
the coupling constant a on the thermalization process. From Eq. (|12j) . we can see that, after the 
mass parameter m is fixed, the event horizon rh is determined by parameters a and b through the 
relation 


(rl + b 


,2 ^ 37-1 _ 


(27) 


with T h = r z + b 2 • When a > 1, there is no bound of b and it can take arbitrary values. When 
a = 1, the above relation reduces to r b + b 2 = 1, so b should be less than 1 to guarantee the 
existence of the horizon. While 0 < a < 1, from Eq. (|14D the maximum value of b is 


bext — ! 'h 


l-3y 


(28) 


The charge parameter q = , ( T> b then takes values in the range q E [0, q ex t\ at this time, with 

Qext = \J 2+a?^ >ex t‘ 

From Eq. (|18D . the ratio of the chemical potential over temperature x is 

47r\/3b 


_h_ = _ 

X ~ T h ^2(2 + a 2 )(r 2 + 6 2 ) 3 / 2 r 


3 ^ [2(37- 1)- 3(2 7 -2)r,]" 1 


(29) 


where we have taken the scale £ = 1. The relation between x and ( q,a ) is plotted in Fig. 1, from 
where we can see that fixing a, 0 < a < 1 and increasing q corresponds to increasing y. When 
q varies from q = 0 to q = q ex ti the value of y varies from 0 to oo. When a > 1, y is no longer 
a monotonically increasing function of q. In fact, as q increases from 0 to oo, y first increases 
and after reaching a maximum value at q m ax it begins to decrease. In this case, the value of y 
is bounded in a finite range. For a = 1, the charge q should be less than \/2 to guarantee the 
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existence of the horizon as we have stated above. From the figure, we can see that in this case, x 
is still a monotonically increasing function of q and approaches a finite constant as q approaches 

V2. 


V/T„ 



FIG. 1 : The chemical potential over the temperature versus q in n = 4 case for various a. For a = 1, q 
should be less than \^2 to guarantee the existence of the horizon. When a > 1, there is a maximum point 
of g/T H at q m ax, which is q max = \/ 6 /( 5 3/8 ) ~ 1.3395 for a = \/2 and q ma x = ~ 1-0745 for a = 2. 


A. The two-point fuction 


On the boundary at time t, we choose two points with spatial separation L and compute the 
correlation function of some operator with large conformal dimension located at the two points. 
Considering the symmetry of the boundary, we can always choose coordinates such that the two 
points lie in the xi-axis with coordinates (t,x i) = (t , — If), (t , ij) respectively. Other coordinates of 
the two points are identical. According to the AdS/CFT duality, in the saddle point approxima¬ 
tion the two-point function can be evaluated by calculating the length of the space-like geodesic 
connecting the two points. The geodesic can be parameterized by only two functions: v(x) and 
z(x), with other coordinates unchanged, where we simply denote x\ = x. Thus, the induced metric 
on the geodesic is 


ds 2 = —N(z)f(v,z)v z — 


N(z) 


2 V 1 + b 2 z 2 


v'z' + 


1 + b 2 z 2 


g(z) dx 2 


(30) 


The length functional of the geodesic is 

r-L/2 


C = 2 


/ P 1/2 dx, 

Jo 


P = — N(z) f (v , z)v' 2 — 


"W *-+ 


1 + b 2 z 2 


2 V 1 + b 2 z 2 


9(z). 


(31) 
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To minimize the length of the geodesic, we need to solve the two equations of motion for v(x) and 
z{x) respectively. The expressions are rather involved and we do not give them explicitly here. 
Considering the symmetry, we have in mind the picture of the geodesic: starting from one of the 
two points on the boundary z = 0, the geodesic extends into the bulk and after reaching the turning 
point (u(0),£(0)) = (u*,z*) it turns back and finally ends at another point on the boundary. At 
the turning point x = 0, the derivatives vanish, that is, t/(0) = z'{ 0) = 0. It seems that one can 
impose the boundary conditions at x = 0 to solve the two equations of motion. However, the two 
equations of motion are singular at x = 0, and in practice, to make our numerical calculations 
available, we would rather impose the boundary conditions at the neighborhood of the turning 
point x = 0, 

v(e) = u* + £>(e 2 ), v'(e) = 0(e 2 ), z(e) = z* + £>(e 2 ), z'(e) = 0(e 2 ), (32) 

where e is a small quantity, with typical order of 10~ 3 . The 0(e 2 ) terms are corrections terms 
and can be derived by solving the two equations of motion around x = 0 order by order. In our 
calculations, we only calculate to order e 2 . The two free parameters w* and z* characterizing the 
turning point of the geodesics are determined by the constraint equations 

v(±L/ 2) = t, z(±L/ 2) = zo, (33) 


where zq is a UV-cutoff to make the length of the geodesic finite and t is the boundary time. 
Furthermore, note that the integrand in the length functional Eq. (|31l) does not depend on x 
explicitly. So if we treat x as a ’’time” variable then the corresponding Hamiltonian is conserved, 
which yields the following equation 


P = 


1 + b z z, 


2 J2 


-g{z*) 


1 + b 2 z 2 


g(z] 


(34) 


By using it, the length functional can be rewritten as 


C = 2 


fL /2 

/ dx 

Jo 


l + b 2 z 2 f 1/2 (1 + b 2 z 2 f , 

— 9 — g{ z *) — 5 — g\ z > 


(35) 


The length of the geodesic £ is a function of the boundary time t. By calculating the length at 
different boundary times, we can get the time evolution of the two-point function. For convenience, 
we define a rescaled length 5C = c ~ Ct} ^ r ™ai an( j s ^ uc iy its time evolution. To make the length 
independent of the cutoff zo, we have subtracted the length at late time in the definition. 6C 
depends on three parameters: the boundary time t , the dilaton-Maxwell coupling constant a and 
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FIG. 2: Time evolution of the two-point function for various values of a and q. Insets in each graph zoom 
in the late time behavior of these curves. The spatial separation of the two points are L = 3. For a = 0, 
q ext = ~ 1.0745; For a = 0.1, q ex t ~ 1.0763; For a = 0.8, q ex t ~ 1.22386; For a = 1, q should be 

smaller than y/2 to guarantee the existence of the horizon. 


the charge q. By varying a and q, we can study the effect of the two parameters on the time 
evolution of the two-point function. 

In Fig. 2, we plot the time evolution of the two-point function for various values of a (0 < a < 1) 
and q. For comparison, we also give the result for the a = 0 case. The spatial separation of the 
two points on the boundary is chosen to be L = 3. From the figure, we can see that that the 
whole process can be divided in the following five successive stages: delay at early times, pre-local- 
equilibrium with quadratic growth in time, post-local-equilibrium with linear growth, memory- 
loss-period, and equilibrium. These features have been observed in various holographic models for 
thermalization and can be considered to be universal for strongly coupled field systems. 

As we increase the charge q, equivalent to increasing the ratio of the chemical potential over 
temperature x °f the dual field system for 0 < a < 1, the initial absolute value of 5C decreases, 
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meaning that the dual field system is initially closer to thermal equilibrium. However, larger 
q makes the delay time longer and the growth of 5C slower, thus the saturation time to reach 
thermal equilibrium increases. This means, for 0 < a < 1, when we fix the energy injected, that 
the saturation time increases as x increases. Also, we should emphasize that the dependence of 
the saturation time on x is weak, as we can see from the figure. This dependence of the saturation 
time on the ratio of the chemical potential over temperature x has been observed in the case with 


standard Maxwell term [23], [2J, |32j, 


and also with non-linear Born-Infeld term 


44]. 


However, when a > 1, the situation becomes complicated, as we show in Fig. 3 with a = 2. In 
this case, the initial absolute value of 5C is no longer a monotonically decreasing function of the 
charge q. In fact, as q increases, the initial absolute value of 5C first decreases and then after q 
crosses some critical value, the initial absolute value of 5C begins to increase. More interesting, 
we can see from the figure that the saturation time is still a monotonically increasing function of 
q. However, from Fig. 1, we can see that for a > 1, x is no longer a monotonically increasing 
function of q, thus the saturation time is no longer a monotonically increasing function of x : When 
q < Qmax i larger x yields longer t sat ; However, when q > q m ax, we have the inverse behavior. We 
can see it more clearly in the right panel of Fig. 3. 




FIG. 3: Left: Time evolution of the two-point function in a = 2 case with L = 3. There is no bound for q 
in this case. q m ax ~ 1.0745 corresponds to the point when reaches its maximum value. There is a large 
overlap between the curves with q = 0.1 and the one with q = 0.8 at early times. Right: Saturation time 
t sa t versus the ratio of chemical potential over temperature \ f° r ct = 2 case. 


By comparing the four graphs with different values of a in Figs. 2 and 3, one can see that the 
effect of the coupling constant a on the saturation time is nearly negligible. We can see this point 
more clearly from the left panel of Fig. 4, where we plot the time evolution of the two-point function 
for various values of a and fixed q. From the figure, we can see that, for fixed q, larger a makes the 
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initial state of the dual field system further away from thermal equlibrium, but will make the delay 
time shorter and the growth of 5C faster thus yielding the saturation time nearly un-affected. The 
saturation time is almost independent of the coupling constant a. As the one-parameter gravity 
action Eq. (JT]) corresponds to a one-parameter dual field theory, this phenomenon indicates the 
universality of the saturation time of this one-parameter strongly coupled field system to reach 
thermal equilibrium under quantum quench. 




L 


FIG. 4: Left: Time evolution of the two-point function for various a with q = 0.8. The spatial separation 
of the two points is L = 3. Right: Saturation time versus the boundary scale. The result is the same for 
various values of a. 

The saturation time depends on the boundary scale L, as we can see from the right panel of 
Fig. 4, where we plot the saturation time versus L. Note that the effect of a on the saturation 
time can be neglected as we stated above, thus this figure is correct for various a. As we can see 
from the figure, the saturation time is linear in L with slope less than unit when L is small, and 
sub-linear when L becomes large. Once again, we can see that the saturation time depends weakly 
on the charge: larger q will have longer saturation time. 

B. The Wilson loop and entanglement entropy 

According to the AdS/CFT duality, the expectation value of the Wilson loop of the boundary 
field theory, in the saddle approximation, is dual to the area of a two-dimensional extremal surface 
in the bulk which anchored on the loop at the boundary. We choose the Wilson loop to be a circle 
and study its time evolution during thermalization. It turns out that the time evolution of this 
observable shows a similar behavior as that of the two-point function, which supports our results 
above. So to avoid redundancy, we will not show the results for the Wilson loop here. 
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Now we turn to use the entanglement entropy to probe the thermalization process. Among 
the three non-local observables, it involves most degrees of freedom of the boundary field system. 
So we expect that by studying its evoluton, we can gain more information of the process. On 
the boundary, we choose a spherical subregion of the spatial space constrained by the relation 
Yt, x 2 -] < R 2 , and calculate its entanglement entropy. According to the conjecture proposed 
by refs. [ZjJ, |&0], the entanglement entropy of a subregion of the boundary spatial space is dual to 
the area of a codimension-two extremal surface anchored on the boundary of the chosen subregion. 
Using spherical coordinates (g, f^) and considering the spherical symmetry, the codimension-two 
extremal surface can be parameterized by only two functions, v(g) and z(g). Thus, the induced 
metric on the surface is 


ds 2 = ~N(z)f(v,z)v' 2 - 


N M „v + 


1 + b 2 z 2 


2 V 1 + b 2 z 2 


g(z) dg 2 + 


1 + b 2 P 


■ g(z)g 2 dnj , ( 36 ) 


and the area functional is 


S = in 

Jo z 


P = -N(z)f(v,z)v' 2 - 


N(z) 

2 V 1 + b 2 z 2 


v’z’ + 


1 + b 2 z 2 


g( z )- 


( 37 ) 


To get the extremal value of the functional, as done in previous two subsections, we need to solve 
the two equations of motion, not shown here. The same boundary conditions as Eq. (13211 are 
imposed. We define a rescaled area 5<S = s ~Stherrnai , w j iere is the volume of the chosen 

subregion, and study its time evolution during thermalization. 

In Fig. 6, we plot the time evolution of the entanglement entropy for various values of a and 
q. The evolution shows a similar behavior as those of the other two non-local observables: for 
fixed a, the saturation time increases as the charge q increases. For 0 < a < 1, this means the 
ratio of the chemical potential over temperature x enhances the saturation time. While for a > 1, 
the situation depends on the charge q: when q < q m ax, the effect of x on the saturation time is 
to enhance, while the effect is opposite when q > q max - Moreover, we see once again that the 
dilaton-Maxwell coupling constant a has almost no influence on the saturation time, although it 
indeed affects the whole process. These results provide further support for the ones obtained from 
the other two observables. 

However, there are some differences comparing to the results of other two non-local observables. 
First, the delay time for the entanglement entropy is the shortest. This is understandable, as 
entanglement entropy, among the three non-local observables involves most degrees of freedom of 
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t 


(c) a = 0.8 


(d) a = 2 


FIG. 5: Time evolution of the entanglement entropy for various a and q. The inset in each graph zooms in 
the late time behavior of each curve. The radius of the entanglement sphere is R = 1.5. 


the field system and thus is the most sensitive to thermalization. Second, now for all cases with 
different a we show, the behavior of the initial absolute value of 5C as q increasing is the same: 
it decreases as q increases. This is different from results of the other two observalbes, where the 
initial absolute value of 5C is not a monotonically decreasing function of q when a > 1. This is 
understandable, as the three observables relate to different amounts of degrees of freedom of the 
boundary field system, and thus may reflect different aspects of the system. As we stated above, 
the entanglement entropy involves more degrees of freedom and thus can can reflect more clearly 
the average state of the field system. Thus it is expected that if all the degrees of freedom of the 
field system are taken into account, the initial state of the field system gets closer to the thermal 
equilibrium as q increases. 
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IV. SUMMARY AND DISCUSSIONS 


We studied the holographic thermalization process in the presence of a chemical potential and 
dilaton field. Three non-local observables are used to probe the process, and they together show 
a rich and consistent picture of the thermalization. The effect of the chemical potential (or the 
charge) and the dilaton field on the thermalization is explored in detail. The saturation time 
increases as the charge q increases, which means that the ratio of the chemical potential over 
temperature x enhances the saturation time when 0 < a < 1. However, when a is larger than 1, 
the situation becomes a little complicated: when q is smaller than some critical value q m axi the 
saturation time increases as the x increases; however, when q > q ma x , the behavior is inverse. This 
is due to the fact that the x is n0 longer a monotonically increasing function of q for a > 1. This 
is different from cases with standard Maxwell term 24, 32] or non-linear Born-Infeld term 44], 
where x is always a monotonically increasing function of q. 

We also investigate the effect of the dilaton-electromagnetic coupling constant a on the thermal¬ 
ization process. This non-trival coupling yields great effects on the availabe background spactimes. 
However, it is remarkable to see that it nearly has no influence on the saturation time. As different 
a in the action defines a different gravity theory and thus corresponds to a different dual held the¬ 
ory, it indicates the universality of the saturation time of this one-parameter held theories to reach 
thermal equilibrium under quantum quench. Although it has no effect on the saturation time, it 
indeed leave imprints on the whole thermalization process. If we hxed the charge q, increasing a 
makes the initial state of the dual held system further away from thermal equlibrium, the delay 
time shorter and the growth of 5C faster. It is interesting to see if an physical explanation on this 
point can be given. 


We have modelled the thermalization process via the naive gravitational collapse of a c 
thin shell. It can also be modelled by the collapsing of a scalar held as done in refs 


larged 


10-14], It is 


interesting to see if this more realistic model can give us more information about the thermalization 
process. We leave it for further investigations. 
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